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I .   NOMENCLATURE 

A.  =  Cross-sectional  area  of  inlet  to  fuel  grain 

Ap  =  Cross-sectional  area  of  fuel  grain 

B  =  Mass  transfer  parameter,  (h^  -  h  ) /h  -  h  ) 

C.  ,C_  =  Coefficients  in  k-e  turbulence  model 

C  =  Specific  heat  at  constant  pressure 

C_  =  Friction  coefficient  without  mass  addition 
fo 

D  =  Inside  diameter  of  fuel  grain 

G  =  Generation  term  in  equations  for  k  and   £  , 


yeff 


M  *  &  •  (W 


/3V    8V  \2< 


G  .  =  Mass  flux  of  air  based  on  bulk  velocity 

air  J 

h  =  Enthalpy,   C   (T  -  T°)  +  m   H  f± 

H  =  Effective  heat  of  combustion,  cal/gm  of  fuel 

i  =  Stoichiometric  ratio  of  oxygen  mass  to  fuel  mass 

k  =  Turbulence  kinetic  energy 

t  =   "Length  scale"  of  eddies  at  near  wall  grid  point  P 

L  =  Length  of  fuel  grain 

m  =  Mass  fraction 

m"  =  Mass  flux  from  fuel  surface 

M  =  Average  molecular  weight 

p  =  Pressure 

P  =  Molecular  Prandtl  number 
r 

r  =  Radial  coordinate 

f  =  Regression  rate  of  fuel  surface 

R  =  Universal  gas  constant 

R  =  Near  wall  grid  Reynolds  number,  p  V  y  /y 


R  =  Internal  radius  of  fuel  port 

1/2 

R  =  Turbulent  Reynolds  number,  It     y  p/y 

S,  =  Source  term  in  differential  equation  for  generalized  dependent 

*  variable  cj> 

T  ■  Temperature 

V  ■  Velocity 

y  =  Distance  from  wall  to  near  wall  grid  point  P 

+  1/2 

y_.  =  Non-dimensionalized  distance  from  wall,  (p/y  )  (t  /p)    y_. 
r                                                           OOP 

z  =  Axial  coordinate 

Y  =  Slip  function  from  Couette  flow  for  boundary  condition  on  k 

£  =  Rate  of  dissipation  of  turbulence  energy 

n.  ■  Combustion  Efficiency  based  on  total  temperature  rise  across 
combustor 

y  =  Molecular  viscosity 

y  rr.  -   Turbulent  viscosity 
err 

p  =  Density 

O   rr  =  Effective  turbulent  Prandtl  or  Schmidt  number 
eff 

T  =  Shear  stress 

cf>  =  Generalized  dependent  variable 

\p  =  Stream  function 

co  =  Vorticity 

Subscripts 

fu  =  Fuel 

h  =  Enthalpy 

N  =  Nitrogen 

o,w  =  Wall 

of  =  Refers  to  the  "property"  mc     -  m  /i 

fu    ox 


ox  =  Oxygen 

P  =  Near  wall  grid  point 

r  =  Radial  component 

T  =  Reference  condition  deep  within  fuel  grain 

z  =  Axial  component 

0  =  Tangential  component 

Superscripts 

o  =  Reference  conditions 

-  =  Average  value 


II.   INTRODUCTION 

Several  characteristics  of  the  solid  fuel  ramjet  indicate  that  it  may  be 
superior  to  other  forms  of  propulsion  for  tactical  weapons  used  at  interme- 
diate ranges  and  high  speed.  Having  no  moving  parts,  the  solid  fuel  ramjet 
is  simple  and  relatively  inexpensive  to  fabricate.  Weight  of  the  system  and 
its  threat  as  a  fire  hazard  are  both  decreased  by  use  of  the  solid  fuel. 

To  be  used  in  a  tactical  situation,  the  solid  fuel  ramjet  has  to  demon- 
strate combustion  stability  and  efficiency  over  the  expected  operating 
envelope  of  altitudes  and  Mach  numbers.   It  must  also  show  performance 
comparable  to  that  of  liquid  fuel  ramjets  and  ducted  rockets. 

Combustion  studies  on  the  solid  fuel  ramjet  have  been  underway  at  United 
Technologies-Chemical  Systems  Division  since  1971.   Initial  work  showed  low 
temperature  rise  combustion  efficiencies.   The  discovery  of  the  rearward 
facing  step  at  the  combustor  entrance  as  a  flameholder  was  a  significant  gain 
in  solid  fuel  ramjet  technology.   Overall  performance,  however,  was  reduced 
by  the  high  stagnation  pressure  loss  produced  by  the  rearward  facing  step. 
Further  work  by  United  Technologies  in  the  field  of  flame  stabilization  in- 
volved various  inlet  designs  including  aerogrids,  distorted  flows,  non- 
circular  inlets  and  vortex  generators.  This  effort  led  to  the  development  of 
inlets  which  minimized  the  required  inlet  step  height  and  decreased  the 
effects  of  inlet  distortion. 

Work  was  also  underway  at  the  Naval  Postgraduate  School  on  the  internal 
ballistics  of  the  solid  fuel  ramjet.   Early  work  showed  that  inlet  turbulence 
and  distortion  may  have  a  significant  effect  on  flame  stability. 

The  solid  fuel  ramjet,  which  uses  air  as  the  oxidizing  agent,  is  similar 
to  the  hybrid  rocket  and  has  two  distinct  combustion  zones  within  the  fuel 
grain  (Fig.  1).  Behind  the  step  is  the  recirculation  zone  where  an  intense 
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mixing  of  reactants  and  products  takes  place.   The  hot  products  ignite  the 
reactant  and  combustion  may  in  the  limiting  case  approach  that  of  a  well- 
stirred  reactor.   This  combustion  region  acts  as  the  flame  initiator  for  the 
combustion  which  takes  place  further  down  the  fuel  grain. 

Downstream  of  the  flow  reattachment  point  a  boundary  layer  develops. 
Here  combustion  is  similar  to  that  of  the  hybrid  rocket.   A  diffusion  flame 
exists  in  the  developing  turbulent  boundary  layer  between  the  fuel  rich  zone 
near  the  wall  and  the  oxygen  rich  central  air  core.  Heat  is  transported  by 
convection  and  radiation  to  the  solid  surface  which  causes  decomposition  of 
the  fuel.   Studies  have  shown  that  the  rate  of  decomposition  of  the  fuel  is 
sensitive  to  the  combustion  pressure,  the  inlet  air  temperature  and  the  mass 
flux  of  the  air.   Finite  rate  kinetics  may  be  of  importance  for  certain 
fuels  and  operating  environments. 

The  combustion  process  in  the  solid  fuel  ramjet  is  thought  to  be  mixing 
limited  at  combustion  pressures  greater  than  10-15  psia  for  all-hydrocarbon 
fuels.   However,  this  remains  to  be  verified  and  other  fuels  have  demon- 
strated behavior  more  characteristic  of  kinetically  controlled  combustion. 

Unburned   gaseous  fuel  escapes  from  under  the  flame  at  the  aft  end  of  the 
fuel  grain  and  results  in  decreased  combustion  efficiency.   This  resulted  in 
the  use  of  a  mixing  chamber  downstream  of  the  fuel  grain.   Temperature  rise 
efficiencies  have  been  found  to  be  a  function  of  the  rate  of  mixing  between 
the  fuel  rich  boundary  layer  and  the  cental  air  core.   Too  rapid  mixing  may 
quench  the  chemical  reactions  and  further  reduce  combustion  efficiency.   Two 
somewhat  distinct  combustion  zones  also  exist  in  the  aft  mixing  section,  the 
recirculation  zone  formed  by  the  fuel  grain  or  aft  orifice  and  the  central 
core  zone0 
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The  use  of  aft-end  mixing  devices  to  improve  combustion  efficiency  is  an 
area  of  recent  interest  in  solid  fuel  ramjet  technology.   As  another  means 
of  promoting  mixing,  the  bypassing  of  a  portion  of  the  inlet  air  around  the 
fuel  grain  and  dumping  it  into  the  aft  mixing  chamber  is  being  re-evaluated. 

A  meteorlogical  sounding  rocket  capable  of  an  altitude  of  200,000  feet  was 

2 
designed  and  built  by  Anderson,  Greenwood  and  Company,  Houston,  Texas  in  1961  , 

The  MET  JET  employed  a  ramjet  using  a  magnesium  and  magnesium-aluminum  alloy 

epoxy-metal  charge.  Eighty-five  percent  of  the  inlet  air  was  bypassed  around 

the  fuel  grain  to  mix  in  an  afterburner  section  with  the  fuel  rich  primary 

flow. 

In  the  bypass  systems  the  flow  rates  into  the  fuel  grain  inlet  and  aft 
mixing  section  are  critical  factors  in  determining  combustion  efficiency. 
Bypass  air  flow  of  too  low  a  percentage  of  the  total  air  mass  flow  or  of  too 
low  momentum  may  have  negligible  effect  on  the  combustion  efficiency.   Bypass 
air  of  too  high  a  flow  rate  or  momentum  may  have  a  negative  effect  on  the 
combustion  process.   The  bypass  air  is  of  appreciably  lower  temperature  than 
the  species  which  exit  from  the  fuel  grain,  and  air  injected  at  too  high  a 
rate  may  cool  the  process  sufficiently  to  affect  the  kinetics  of  the  reaction. 
Furthermore,  the  combustion  process  in  the  aft  mixing  chamber  may  also  be  a 
function  of  the  axial  and  radial  positions  and  the  angular  orientations  of 
the  aft  dumps. 

This  experimental  investigation  considered  the  effects  of  bypass  con- 
figuration and  operating  conditions  on  the  combustion  behavior  of  a  solid  fuel 
ramjet  which  utilized  polymethylmethacrylate  as  the  fuel. 

In  order  to  guide  development  efforts  and  to  provide  needed  input  to 
systems  analysis  studies,  combustion  models  are  needed  which  can  predict  the 
fuel  regression  rate,  flammability  limits  and  combustion  efficiency  as  a 


function  of  hardware  design  and  operating  environment. 

Modeling  efforts  for  the  combustion  process  in  solid  fuel  ramjets  haye 

3 
proceeded  along  two  directions.   Dunlap  has  developed  a  semi-empirical 

regression  rate  model  which  is  based  upon  reacting  turbulent  boundary  layer 

theory  and  experimentally  obtained  wall  heat  transfer  rates  in  non-reacting 

1  4 
flows.  Work  at  the  Naval  Postgraduate  School  *   has  been  directed  toward 

adaptation  of  the  Spalding,  et  al  model  '   for  heat  and  mass  transfer  in 

recirculating  flows  to  the  SFRJ  geometry. 

Experimental  data  3  have  shown  that  the  recirculation  zone  is  fuel  rich 
and  that  the  primary  combustion  region  (or  flame  zone)  spreads  out  from  along 
the  shear  layer  between  the  recirculation  zone  and  the  inlet  flow.   The 
flame  has  been  observed  to  be  quite  broad  near  the  flow  reattachment  posi- 
tion and  subsequently  develops  into  a  turbulent  diffusion  flame  within  the 
developing  boundary  layer  downstream  of  reattachment. 

The  initial  modeling  work  at  the  Naval  Postgraduate  School  considered 
the  effects  of  finite  rate  kinetics  and  combustor  geometry  on  the  flame 
pattern  and  internal  flow  field. 

Several  weaknesses  were  evident  in  the  earlier  work.   Fuel  regression 
rates  were  not  calculated.   Also,  the  flame  pattern  was  not  in  agreement 
with  experimental  data,  i.e.,  the  recirculation  region  was  predicted  to  be 
oxidizer  rich  and  the  flame  within  the  boundary  layer  region  spread  too 
rapidly  toward  the  centerline  of  the  motor.   In  addition,  the  model  had 
not  been  adequately  checked  for  accuracy  against  experimental  data  for 
velocity,  pressure  and  turbulence  intensity. 

The  present  investigation  was  conducted  to  improve  the  qualitative 
accuracy  of  the  model,  to  incorporate  prediction  of  the  fuel  regression 
rate,  and  to  check  the  validity  of  the  model  against  experimental  data. 
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The  majority  of  the  analytical  studies  were  conducted  for  polymethyl- 
methacrylate (PMM)  fuel.   PMM-air  combustion  produces  minimal  radiative  heat 
transfer  to  the  fuel  grain  for  high  mass  flow  rates  through  the  port.   The 
latter  condition  is  characteristic  of  non-bypass  systems.   The  model 
developed  for  this  system  did  not  include  radiative  heat  transfer  or  finite 
rate  kinetics. 

The  current  investigation  also  considered  the  application  of  the  model 
to  the  combustion  behavior  of  all-hydrocarbon  fuels  which  exhibit  significant 
gas-phase  radiation.   In  particular,  model  predictions  were  compared  to  the 
data  obtained  by  Schadow  . 

III.   EXPERIMENTAL  INVESTIGATION 
A.   Equipment  and  Procedures 

A  schematic  of  the  solid  fuel  ramjet  is  shown  in  Fig.  2.   The  apparatus 
employed  a  fuel  port  to  dump  inlet  area  ratio  of  9.0  and  a  fuel  port  to 
nozzle  throat  area  ratio  of  4.0.   The  exit  area  of  the  grain  (for  most  tests 
conducted)  was  held  fixed  at  the  initial  port  area  by  using  a  thin  orifice 
plate.   The  aft  mixing  chamber  had  a  length  to  diameter  ratio  (L/D)  of  2.93. 
The  aft  mixing  chamber  consisted  of  three  interchangeable  sections  such  that 
the  axial  location  of  the  bypass  dump  could  be  varied.   One  section  contained 
four  inlets  for  the  bypass  air.   The  inlet  diameters  were  2.70  cm  with  plugs 
available  to  reduce  the  diameters  to  2.065,  1.908,  1.458,  or  0.635  cm,  or  to 
seal  the  inlets  off  entirely.  An  additional  section  was  fabricated  which  pro- 
vided bypass  dumps  oriented  tangentially  to  the  mixing  chamber  wall  and 
perpendicular  to  its  centerline.   In  most  tests  conducted  the  bypass  dump  was 
located  in  the  forward  position,  within  the  recirculation  zone  of  the  aft 
mixing  chamber. 
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All  test  firings  were  performed  in  the  jet  engine  test  cell  at  the  Naval 
Postgraduate  School.   Data  for  calculating  temperature  rise  efficiencies  based 
on  combustion  pressure  were  obtained  while  varying  primary  and  bypass  air  flow 
rates  and  bypass  dump  geometry.   Nominal  air  flow  rates  were  0.1  and  0.2  lbm/ 
sec  for  non-bypass  runs  and  0.2  lbm/sec  for  runs  using  bypass.   A  summary  of 
the  test  conditions  are  presented  in  Table  I. 

Table  I  Summary  of  Experimental  Conditions 

Air  Inlet  Temperature  (°K) :   292,  389 

Bypass  Dump  Location:   forward  and  aft  of  flow  reattachment 

Bypass  (%  Port/%  Bypass):   100/0,  65/35,  50/50.  35/65 

Momentum  Ratio  (Dump /Port) 

for  50/50  Bypass:   0.3,  0.5,  0.6,  1.0,  5.3 

Dump  Configuration:  none,  2  at  180°,  2  at  90°,  4  at  90°,  2  at  180°  with  swirl 

The  majority  of  the  test  firings  were  made  using  a  nominal  inlet  air 
temperature  of  70°F.  For  those  using  the  air  heater  (non-vitiated  air),  two 
to  three  hours  were  allowed  for  the  temperature  to  stabilize  at  the  ramjet 
inlet  and  in  the  bypass  air  line. 

The  weight,  length  and  inside  diameter  of  the  fuel  grains  were  measured 
prior  to  being  mounted  in  the  motor.   The  ignition  sequence  normally  lasted 
for  six  seconds.   An  average  of  three  seconds  was  required  for  the  ignition 
flame  to  propagate  from  the  head-end  assembly  into  the  fuel  grain.   After 
five  seconds,  the  air  was  directed  through  the  ramjet  motor.   Ignition  was 
continued  into  the  first  second  of  the  run  to  insure  that  combustion  would 
be  sustained.   Total  time  of  the  ignition  flame  in  contact  with  the  fuel 
grain  amounted  to  approximately  three  seconds.   Two  tests  were  made  using 
only  the  ignition  system  in  order  to  determine  the  rate  of  consumption  of 
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the  PMM  grains  during  the  oxygen-methane  ignition.   These  data  were  used  to 
correct  the  initial  weight  of  the  fuel  used  in  the  efficiency  calculation. 

Combustion  normally  lasted  for  forty-five  seconds.   The  motor  was  ex- 
tinguished at  the  end  of  each  run  by  simultaneously  venting  the  air  to  the 
atmosphere  and  actuating  the  nitrogen  purge  system.  Low  pressure  air  was 
then  blown  through  the  motor  for  cooling. 

Thirty-three  hot  firing  PMM  tests  were  conducted.   Combustion  efficien- 
cies were  based  on  the  total  temperature  rise  from  just  upstream  of  the  rear- 
ward facing  step  at  the  fuel  port  inlet  to  the  nozzle  inlet.   Inlet  total 
temperature  was  derived  from  the  measured  inlet  static  temperature.  Actual 
combustion  temperature  was  derived  from  the  measured  static  pressure  in  the 
aft  mixer o 

Regression  rate  was  calculated  based  on  weight  loss  of  the  fuel.   The 
Inside  diameter  of  the  aft  end  of  the  fuel  grain  was  also  measured  before  and 

o 

after  each  firing.   However,  as  observed  by  Boaz  and  Netzer  ,  it  was  found 
that  weight  loss  gave  a  more  consistent  value  of  regression  rate  than  the 
method  based  on  aft  end  diameter  change,  due  to  the  non-uniform  regression 
along  the  length  of  the  grain. 

The  calculated  temperature  rise  efficiencies  based  on  pressure  are  prone 
to  error,  since  uncertainty  in  efficiency  varies  as  the  square  of  the  uncer- 
tainty in  the  pressure  measurement.   Errors  as  large  as  7%  were  possible. 
However,  the  uncertainty  in  percent  theoretical  C*  was  approximately  2%  and 
the  temperature  rise  efficiencies  followed  the  C*   efficiencies  within  1.5%. 
B.   Results  and  Discussion 

(a)   Regression  Rate 

Several  runs  were  made  without  using  bypass  air.   These  showed  the  depen- 
dence of  the  regression  rate,   f  (in/sec)   on  combustion  pressure,   P  (psia) 
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2 
and  average  mass  flux  of  air,   G  .   (lbm/ia  -sec)  to  be: 

29 
f  =  0.0043P*  yG    .    .38  (1) 

air  v  ' 

A  plot  of  the  regression  rate  versus  this  empirical  regression  rate  equation 
for  the  0.2  and  0.1  lbm/sec  non-bypass  runs  made  with  cold  inlet  air  is 

shown  in  Fig.  3. 

g 
In  their  work  with  PMM,  Boaz  and  Netzer  had  shown  a  regression  rate 

dependence  of  the  form: 

r  =  CP*51T  .  *34G  .  *41  (2) 

air    air  ' 

Insufficient  runs  were  made  in  this  investigation  at  other  than  ambient 
inlet  temperature  to  calculate  a  temperature  dependence  for  the  regression 
rate.   The  above  equations  for  regression  rate  agree  closely  for  the  depen- 
dence on  G.   The  dependence  on  pressure  was  significantly  less  in  this  study. 
In  the  work  of  Boaz  and  Netzer,  chamber  pressure  was  varied  intention- 
ally from  37  to  108  psia  by  varying  the  throat  size.   In  this  study  pressure 
varied  between  33  and  63  psia  only  as  a  result  of  varying  combustion  effi- 
ciency and  G.   The  dependence  of  regression  rate  on  G  results  from  convective 
heat  transfer  to  the  fuel  surface  whereas  the  dependence  on  the  pressure  can 
result  from  radiation  and/or  finite  rate  kinetics. 

With  the  application  of  bypass  air,  the  dependence  of  regression  rate  on 
pressure  and  air  flux  was  altered.   In  Fig.  4  are  plotted  the  0.2  lbm/sec 
non-bypass  cases  and  the  bypass  runs  made  with  0.2  lbm/sec  total  air  flow  rate 
at  the  nominal  test  condition.   The  nominal  test  condition  was  a  bypass  air 
flow  rate  of  0.1  lbm/sec  and  two  dumps  with  diameters  of  2.065  cm.   This 
corresponded  to  a  dump  momentum  to  fuel  grain  port  momentum  ratio  of  approxi- 
mately 0.5.  A  slightly  stronger  dependence  on  pressure  was  shown  while 
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regression  rate  indicated  very  little  or  no  dependence  on  air  flux  for  the 
values  tested.  Here  regression  rate  took  the  form: 

r  =  0.00116P,42G  .  •003  (3) 

air 

In  the  bypass  situation,  the  mass  flux  through  the  grain  is  low  but  the 
pressure  is  maintained  high  due  to  the  total  mass  flux  through  the  nozzle 
throat.  However,  correcting  the  regression  rate  for  the  increased  pressure 
per  equation  (2)  does  not  result  in  regression  rates  as  high  as  the  experi- 
mental data.   In  addition,  regression  rates  based  on  weight  and  diameter 
agreed,  indicating  that  the  change  did  not  result  from  different  combustion 
behavior  within  the  aft  mixing  chamber.   These  conditions  of  low  G  and  high 
P  minimize  the  convective  heat  flux.  At  lower  G  the  regression  rate 
increases  relative  to  the  air  flux,  i.e.,  more  unburned  fuel  exists  under- 
neath the  diffusion  flame  within  the  thicker  boundary  layer;  thus  more  gas 
with  radiative  properties  is  present. 

(b)   Combustion  Efficiency 

In  Fig.  5  are  plotted  the  cases  from  both  Figs.  3  and  4.  As  can  be 
seen  from  this  figure  the  use  of  bypass  air  flow  in  a  solid  fuel  ramjet  using 
polymethylmethacrylate  as  a  fuel  has  the  effect  of  reducing  combustion 
efficiency.   For  the  non-bypass  runs,  a  decrease  in  air  flow  also  causes  a 
slight  decrease  in  combustion  efficiency.  While  maintaining  the  same  air 
flux  through  the  grain,  injecting  bypass  air  into  the  mixing  section  brings 
a  further  decrease  in  performance.   Decreasing  air  flux  through  the  fuel 
grain  while  maintaining  the  same  total  flow  rate  in  the  bypass  case  also 
brings  about  a  reduction  in  combustion  efficiency.   Decreasing  the  air  flux 
through  the  fuel  grain  increases  the  percentage  of  fuel  that  must  be  burned 
in  the  aft  mixing  chamber.  Apparently,  the  mixing  chamber  was  not  of 
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sufficient  length  and/or  the  bypass  air  caused  quenching  of  the  mixing 
chamber  combustion. 

To  further  study  the  decrease  in  performance  found  when  using  bypass, 
other  test  conditions  were  considered.   Two  runs  were  made  with  the  dumps 
located  behind  the  aft  mixing  chamber  reattachment  point.   There  was  only  a 
slight  decrease  in  combustion  efficiency  (81%)  over  that  of  identical  runs 
in  the  forward  position  (83%) . 

Three  runs  were  made  varying  the  momentum  of  the  individual  dumps  from 
the  nominal  test  condition.   Small  changes  in  momentum  had  no  effect.  Vary- 
ing the  momentum  ratio  from  0.3  to  1.0  did  not  significantly  change  the  effi- 
ciency. However,  when  a  significantly  large  increase  in  momentum  was 
effected  (5.3),  the  combustion  performance  showed  a  noticeable  decrease  (73%). 

Only  limited  data  were  taken  for  other  than  two  bypass  dumps  at  180 
degrees.   However,  the  data  showed  that  with  low  dump  momentum  two  or  four 
dumps  spaced  at  90  degrees  reduced  the  combustion  efficiency.   Four  dumps  at 
ninety  degrees  with  high  momentum  decreased  the  combustion  efficiency  slightly, 

The  last  variation  in  geometry  studied  was  the  use  of  swirl  in  the  bypass 
dump  process.   This  was  accomplished  by  injecting  the  bypass  air  with  a  tan- 
gential velocity  component.   The  solid  fuel  ramjet  motor  was  also  run  without 
fuel  grain  ignition,  both  with  and  without  swirl,  to  determine  the  effect  of 
the  induced  vorticity  on  the  effective  exit  nozzle  diameter.   Decreasing  the 
effective  nozzle  diameter  would  increase  the  combustion  chamber  pressure  and 
give  false  indications  of  higher  efficiency.   The  swirl  was  found  to  not 
affect  the  effective  throat  diameter.   Two  test  firing  runs  were  made  with 
swirl,  giving  an  average  combustion  efficiency  of  86%,  the  highest  perfor- 
mance of  any  bypass  run. 
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Every  form  of  bypass  used  caused  a  decrease  in  combustion  efficiency. 
This  indicates  that  the  decrease  in  performance  is  due  to  the  effect  of  the 
bypass  air  flow  on  the  kinetics  of  the  combustion  process  within  the  aft 
mixing  chamber.  Essentially  equal  decreases  in  performance  were  obtained 
with  the  dump  air  injected  both  behind  and  in  front  of  the  reattachment  point. 
Apparently  a  significant  portion  of  the  total  combustion  process  takes  place 
downstream  of  the  reattachment  point  in  the  aft  mixing  chamber.   This  agress 
with  the  temperature  data  presented  by  Schadow  for  an  all-hydrocarbon  fuel. 
In  the  case  of  PMM,  the  light  weight  unburned  hydrocarbons  that  enter  the  aft 
mixing  chamber  apparently  burn  most  completely  when  allowed  to  react  slowly 
with  the  available  oxygen  in  the  hot  flow  from  the  core  of  the  fuel  grain. 

The  possibility  of  the  temperature  of  the  air  entering  the  aft  mixer 
causing  reduced  combustion  efficiency  was  considered.  A  run  was  made  in  which 
the  temperature  entering  both  the  grain  and  aft  mixer  was  increased  from 
approximately  292°K  (520°R)  to  389°K  (750°R).   The  combustion  efficiency  did 
not  change  significantly,  indicating  that  for  PMM-air  combustion  the  dump  air 
temperature  is  probably  not  as  important  as  the  quantity  and  temperature  of 
the  unburned  fuel  and  the  mixing  rate. 

In  the  region  in  front  of  the  reattachment  point,  the  less  that  was  done 
to  disturb  the  flow  resulted  in  better  performance.   In  the  case  of  swirl 
where  the  bypass  flow  remained  close  to  the  wall,  bypass  efficiency  was  maxi- 
mized. When  the  bypass  momentum  was  increased  to  where  the  flow  disturbed  the 
fuel  rich  layer  between  the  recirculation  region  and  the  air  rich  central  core, 
performance  decreased.   This  again  indicates  that  a  major  portion  of  the  com- 
bustion process  takes  place  along  this  fuel  rich  layer  or  that  the  process 
downstream  is  highly  dependent  on  the  high  temperatures  in  this  layer.   It 
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also  indicates  that  the  combustion  mechanisms  around  this  layer  are  more 
important  than  those  which  occur  within  the  recirculation  region.   These  re- 
sults imply  that  different  optimum  bypass  dump  configurations  and  momentum 
should  be  expected  for  different  fuel  systems. 

It  is  known  that  the  use  of  bypass  improves  performance  in  solid  fuel 
ramjets  using  all-hydrocarbon  fuels.   In  the  case  of  oxygen-containing  fuels 
and  for  fuels  which  decompose  into  monomers  or  small  hydrocarbon  molecules 
(such  as  PMM) ,  results  indicate  that  the  fuel  burns  most  efficiently  without 
bypass. 

The  use  of  bypass  systems  has  meant  an  increase  in  weight,  cost  and  com- 
plexity of  the  solid  fuel  ramjet.   In  addition,  they  may  introduce  combustor- 
feed  system  coupling.   The  use  of  a  fuel  which  has  sufficient  density  impulse, 
regression  rate  and  flammability  limits  to  minimize  inlet  total  pressure 
losses  has  led  to  the  use  of  all-hydrocarbon  fuels.   Although  PMM  does  not 
meet  the  criteria  for  a  good  fuel,  the  results  of  this  study  indicate  that 
future  fuel  studies  may  be  fruitful  if  directed  toward  ones  which  contain  low 
percentages  of  oxidizer  and/or  substances  which  unzip  the  hydrocarbon  chain. 

(c)  Combustion  Pressure  Oscillations 

In  the  non-bypass  test  runs,  the  inlet  and  combustion  pressures  exhibited 
a  steady,  small  amplitude  (approximately  2%  of  chamber  pressure)  oscillation 
of  approximately  150  Hz.   In  the  bypass  runs  at  nominal  conditions,  the  same 
oscillation  appeared  though  of  considerably  higher  amplitude  (approximately 
30%  of  chamber  pressure). 

In  the  case  of  low  or  very  high  aft  dump  momentum,  a  second  oscillation 
appeared  along  with  that  previously  mentioned.   This  oscillation  was  of  very 
low  frequency  (1  Hz.)  and  large  amplitude  (approximately  20%  of  chamber  pres- 
sure) and  may  be  connected  with  behavior  within  the  mixing  chamber  recirculation 

18 


zone.   The  bypass-generated  very  low  frequency  pressure  oscillations  are  in 
the  range  of  frequencies  which  could  be  coupled  to  the  fuel  regression  rate. 
Other  fuel/feed  systems  may  exhibit  similar  or  different  oscillatory  behaviors 
depending  upon  the  feed  system  design  and  the  kinetics  of  the  combustion 
process  in  the  aft  mixing  chamber. 
C.   Conclusions  and  Current  Work 

The  use  of  bypass  air  flow  in  a  solid  fuel  ramjet  which  utilizes  PMM  as 
the  fuel  decreases  the  combustion  efficiency.   Air  flow  injected  into  the  aft 
mixing  chamber  has  a  more  pronounced  effect  on  the  combustion  process  when  a 
high  enough  momentum  is  provided  for  the  bypass  air  to  reach  the  fuel  rich 
shear  layer  trailing  from  the  port  of  the  fuel  grain. 

A  significant  amount  of  combustion  occurs  downstream  of  the  reattachment 
point  in  the  aft  mixing  chamber.  Bypass  air  injected  into  the  aft  mixing 
chamber  also  has  an  effect  on  the  regression  rate  upstream  in  the  fuel  grain. 
For  PMM  fuel  burned  at  high  pressures  and  low  air  mass  flux,  the  principal 
mechanism  for  wall  heat  flux  became  radiation,  and  resulted  in  the  regression 
rate  becoming  insensitive  to  the  air  flux  through  the  fuel  port. 

With  the  use  of  bypass  systems  it  is  possible  to  set  up  combustor-f eed 
system  type  oscillations  and  instabilities  dependent  upon  the  effects  of  the 
bypass  flow  on  the  aft  mixing  chamber  combustion  process  and  on  the  fuel 
regression  rate. 

Current  experimental  work  is  being  directed  at  the  effects  of  bypass  on 
the  combustion  behavior  of  all-hydrocarbon  fuels.   In  addition,  the  effects 
of  design  and  test  conditions  on  combustion  pressure  oscillations  are  being 
investigated. 


19 


IV.   ANALYTICAL  INVESTIGATION 
A.   Model  Overview 

The  basic  assumptions  were  that  the  flow  was  steady,  recirculating, 
two-dimensional  and  subsonic.   In  addition  kinetic  heating  and  vorticity 
sources  resulting  from  spatial  gradients  in  the  effective  viscosity  were 
neglected.   The  turbulent  Lewis  number  was  taken  to  be  unity.   The  trans- 
formed variables  of  vorticity  (w/r)   and  stream  function  (¥)  were  used 
together  with  enthalpy   (h)  ,  and  species  mass  fractions  (oxidizer,  fuel, 
etc.).   A  modified  Jones-Launder  •  '    two  equation  turbulence  model  was 
used  to  calculate  the  viscosity  throughout  the  flow  field.   Thus,  two  addi- 
tional variables  must  be  considered,  the  turbulence  kinetic  energy  (k)   and 
the  turbulence  energy  dissipation  rate  (e)  .  The  modified  Jones-Launder 
model  incorporates  five  empirical  constants  and  the  effective  viscosity  is 
calculated  using  the  following  expression: 

Ueff  =  0-°9pk2/e  (4) 

In  the  present  study  the  combustion  was  considered  to  be  mixing 
limited  with  a  one-step  simple  chemical  reaction.   Thus,  1  gm  fuel  +  i 
gm  oxidizer  ■*■     (1+i)  gm  of  products  and  fuel  and  oxidizer  could  not  exist 
together.   Species  considered  were  oxygen,  nitrogen,  fuel  and  products. 
The  present  form  of  the  model  neglected  radiation  and  density  was  calculated 
using  p  =  pM/RT  . 

With  these  assumptions  the  governing  equations  for  axisymmetric  flow 
can  be  cast  in  the  standard  elliptic  format  '  ' 

r3      ,.    34\        3      ,.    W,_        3      ru        HC^\        3.     rv  V1  B  „     ,« 

aA    ["5-  (<$>  -K-)    -  T~  (4>  -st-)]    "  T"   [b.r  — r— ]    -  ■%-  [b.r  — s^ — ]   +  r   S  ,    =   0      (5) 

(f)      3z     v   3r  3r     Y   3z  3z       tp  3z      J        3r        <J>  3r  <j> 

convection  diffusion  source 

terms  terms  term 
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where  the  parameters  are  presented  in  Table  II.  Additional  discussion  of 
these  equations  and  utilization  of  various  turbulence  models  can  be  found 
in  Ref.  11  and  12. 


Table  II.   Equation  Parameters 
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B.   Boundary  Conditions  And  Solution  Procedure 

The  inlet  flow  was  considered  to  be  "plug  flow"  although  other  inlet 
conditions  (distortion,  etc.)  are  easily  incorporated.   The  inlet  turbu- 
lence intensity  was  taken  to  be  5.2%  in  order  to  closely  match  experimental 
measurements.   Axial  gradients  for  all  variables  were  considered  zero  at 

the  flow  exit  and  were  taken  as  zero  for  all  "conserved  properties" 

m 
(mf t-   ,  m^  ,  h)   at  the  step-face.   Radial  gradients  for  all  varia- 
bles except  ¥  were  taken  as  zero  at  the  axis  of  symmetry.   The  major 
area  of  concern  in  the  specification  of  boundary  conditions  for  the  SFRJ  is 
along  the  fuel  surface.  When  calculations  are  made  with  the  model  for  an 
adiabatic  or  isothermal  wall  with  no  mass  addition,  inaccuracies  in  the 
turbulence  model  are  somewhat  masked „   The  fuel  mass  addition  at  the  wall 
and  the  subsequent  combustion  greatly  amplify  the  need  for  accurate  effec- 
tive viscosity  calculations  near  the  wall  if  observed  experimental  behavior 
is  to  be  approximated. 

All  "conserved"  properties  when  in  dimensionless  form,  had  identical 
governing  equations  and  boundary  conditions.   Thus,  only  one  equation  was  re- 
quired to  be  solved  for  all  "conserved"  properties.   The  resulting  set  of  five 
equations  (— ,  \pt   k,  e,  m   -  m  /i)  were  cast  in  finite  difference  form  and 
solved  using  the  Gauss-Seidel  method  with  upwind  differences  and  relaxation. 
Under  relaxation  was  required:  0.5  for  —  and  ijj  ,  0.3  for  k  and  £  ,  and 
0o9  for  mf  -  m  /i  .   A  non-uniformily  spaced  17  x  25  rectangular  grid  system 
was  employed.   Solution  was  performed  on  an  IBM  360-67  using  Fortran  H  and 
required  approximately  thirty-five  minutes  to  complete  1000  iterations. 

The  fuel  surface  was  considered  to  have  a  uniform  and  constant  temperature 
and  the  mass  fraction  of  oxidizer  at  the  surface  was  considered  to  be  zero. 
In  order  to  obtain  a  fuel  rich  recirculation  zone,  a  boundary  layer  flame 
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pattern  in  agreement  with  experiment,  and  accurate  conservation  of  species, 
several  modifications  had  to  be  made  to  the  original  model  '  '  .   For  simpli- 
city the  boundary  layer  was  considered  to  consist  of  two  parts,  a  laminar 
sublayer  and  a  turbulent  layer.   The  border  between  the  two  regions  was 
considered  to  occur  where  y  =  24  .  At  each  near  wall  grid  point  y   was 
calculated. 

If  y+  >  24:   Cr  II   =  0.02846/R  °'25  ($) 

p        fo  e  v  ' 

13 
The  velocity  profile  was  corrected  for  blowing   and  u  „_  was  calculated 

ef  f 

using  Eq.  (4) . 

If  y+  <  24:   C,  /2  =  2/R  .  (7) 

'p  —      fo       e  v  ' 

The  laminar  viscosity  was  calculated  using 

V   -  PQ  (T/To)0'5  (8) 

In  the  actual  calculation  procedure,  the  grid  was  selected  to  insure 
that  y  <  24  .   The  value  of  0.02846  is  required  for  matching  of  the 
two-zone  boundary  layer. 

The  boundary  condition  for  vorticity  was  calculated  from 

T  2 

(<£)   =  _° —  where  x  =  (C£  IT)    pV  In  (1  +  B) /B        (9) 
r  o   u  r  o     fo 

o  w 

in  which  the  wall  shear  stress  is  corrected  for  blowing  by  using  the  simple 
Couette  flow  approximation  ln(l  +  B) /B  .   The  control  volumes  used  at  the 
fuel  surface  for  the  micro-integration  are  shown  in  Fig.  6.   The  control 
volume  for  all  variables  but  k  and  £  were  taken  to  the  wall.   For  k 
and  e  the  control  volumes  were  taken  to  the  half -grid  point.   Gradients 
at  the  wall  for  "conserved"  properties  were  taken  to  be  of  the  form 
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(|£)   -  {(h  -  h  )/y  }  In  (1  +  B)/B  (10) 

dy  o      p    op 


The  wall  mass  flux  provided  the  boundary  condition  for  ¥  and  was 
calculated  from 

;..  =  £  0  /(h  .  v  =  (Ja,  iniia!  (11) 

rr    '  *r     'p 

The  boundary  condition  on  k  was  calculated  from 

k  =  y  k  (12) 

o   '   p  ' 

Y  -  -  1.0  ,  R  <  19.8 

=  -  0.39,  R  >  19.8 
t 

The  boundary  condition  on  z     was  considered  to  occur  at  the  near 

wall  grid  point 

S   -  kp3/2  '\  W 

where  i     =0.2y   ,  y  +  <  24  (14) 

P       P  '   P 

I     =  0.4  y  ,  y  +  >  24  (15) 

P       P   'P 

The  boundary  conditions  for  —  ,  k  and  z     on  the  step  face  were 

calculated  in  a  similar  manner  except  that  there  was  no  wall  mass  addition, 

none  of  the  micro-integration  control  volumes  were  extended  to  the  wall,  and 

I     =   2.44  y   . 

P         P 

14 
Launder  and  Spalding   have  discussed  the  strengths  and  weaknesses  of 

the  k-£  model  using  both  the  wall  function  and  low-Reynolds-number  modeling 

methods.   They  indicate  that  modification  of  the  boundary  conditions  on  E 

were  sometimes  found  necessary  on  the  step  face  in  order  to  obtain  results  in 

agreement  with  experiment.   In  this  investigation  wall- functions  were  employed 
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for  the  k-e  model.   Since  the  mass  addition  rates  were  low  (B  <  1.5)  the 

five  constants  in  the  model  *  '    were  not  modified.  When  y   <  24  along 

P 

the  fuel  grain  the  effective  viscosities  predicted  using  the  k-e  model  were 
not  used  at  the  near  wall  grid  points.   In  this  case,  the  boundary  conditions 
given  by  Eqns.  (13),  (14)  and  (15)  were  required  to  provide  realistic  values  of 
u  ff  and  fuel  mass  fraction  at  the  second  grid  point  from  the  wall.  An  alter- 
nate procedure  could  have  been  to  use  the  low-Reynolds-number  modeling 
methods  '  .   In  recent  work  with  smaller  inlet  step  heights  it  has  been  found 
necessary  to  employ  Eqns.  (13),  (14),  and  (15)  for  e   along  the  step  face. 
Much  additional  work  is  apparently  required  to  better  define  the  near-wall 
region. 

The  finite-difference  equations  '  require  the  use  of  space  averaged 
values  for  the  source  terms  presented  in  Table  II.  At  the  near  wall  grid 
points  large  gradients  exist  if  y  >  24  .  As  mentioned  above,  this  was 
prevented  at  the  fuel  surface  by  appropriate  selection  of  the  radial  grid 
spacing.  Along  the  step  face,  however,  integrated  average  values  had  to  be 
calculated  for  —  and  all  velocity  gradients.   Velocity  profiles  along  the 

step  face  for  y   >  24  were  assumed  to  obey  the  l/7th  power  law. 
P 

In  order  to  obtain  realistic  solutions  the  stream  function  at  the  grid 
point  downstream  of  the  step  inlet  was  required  to  be  held  fixed  at  the 
value  which  existed  at  the  step  corner. 
C.   Results  and  Discussion 

Divergence  would  occur  in  the  iteration  unless  the  initial  conditions 
specified  for  co/r  and  i[»  were  in  close  agreement  (for  example,  l/7th 
power  law  distributions).   It  was  also  necessary  to  prescribe  an  initial 
length  scale  distribution.   The  length  scale  distribution  was  used  to  cal- 
culate initial  values  of  z     from  k  using  Eqn  (13) .   The  £     field  was 
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held  constant  for  an  initial  35  iterations  in  order  to  promote  convergence. 
Specification  of  the  initial  length  scale  distribution  required  trial  and 
and  error  methods  for  other  than  the  simple  rearward  facing  step  geometry. 
Variations  in  gas  density  were  also  prevented  from  varying  for  the  first  50 
iterations. 

At  each  axial  location  species  and  energy  conservation  calculations 
were  made.   False  diffusion  errors  were  larger  in  the  recirculation  zone 
where  the  grid  lines  were  at  large  angles  to  the  streamlines.   Small  errors 
in  the  solution  of  the  "conserved"  property  equation  presented  some  difficul- 
ties since  mf   ,  m   ,  m    and  m.,   were  all  calculated  from  the  one 
fu  '   ox    pr        N2 

solution.   Calculation  of  fuel  conservation  was  especially  difficult  since 
fuel  concentration  gradients  near  the  wall  were  very  large.   In  general, 
energy  and  species  (other  than  fuel)  were  conserved  to  within  five  percent. 
With  these  accuracies,  fractional  changes  between  iterations  were  a  maximum 
of  0.3  percent.   This  required  between  750  and  1000  iterations.   The  "conserved' 
properties  and  stream  function  converged  more  rapidly  than  did  —  ,  k  and  e    , 
For  the  higher  inlet  velocities,  the  relaxation  parameter  for  —  had  to  be 
reduced  to  values  between  0.3  and  0.4. 

The  location  of  the  "reattachment  point"  was  found  to  be  somewhat  sen- 
sitive to  the  distance  from  the  step  face  to  the  first  axial  grid  line 
This  resulted  from  specification  of  the  stream  function  as  unchanging  as 
discussed  above.   However,  once  a  particular  grid  was  chosen  which  yielded 
a  reattachment  point  in  agreement  with  experiment  for  one  step  height  and 
inlet  velocity,  good  results  were  obtained  for  other  step  heights  and  inlet 
conditions.   The  reattachment  point  was  found  to  be  insensitive  to  inlet 
velocity  and  was  found  to  move  slightly  upstream  with  low  rates  of  wall  mass 
addition.   Both  of  these  results  are  in  good  agreement  with  experimental 

data 
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A  true  reattachment  point  does  not  exist  when  mass  addition  through 
the  wall  is  present.  However,  a  pseudo-reattachment  point  was  calculated 
by  assuming  it  to  occur  where  the  average  velocity  was  zero  between  the  wall 
and  the  near-wall  node.  A  comparison  of  this  "reattachment  point"  with 
chemical  reaction  is  made  with  experimental  data  from  non-reacting  flows  in 
Fig.  7.   It  should  be  noted  that  the  predicted  flow  reattachment  occurred 
downstream  of  the  Krall-Sparrow1  maximum  heat  transfer  point  but  upstream 
of  the  Phaneuf-Netzer1^  data  with  10%  wall  mass  addition. 

Typical  predicted  streamline  patterns  and  temperature  distributions 
are  presented  in  Figure  8  for  plexiglas  fuel  and  air.  The  dividing  stream- 
line is  also  sketched  in  the  temperature  field.   The  recirculation  region 
was  calculated  to  be  fuel  rich  and  the  flame  initiated  along  the  shear  layer 
and  spread  out  as  it  approached  the  reattachment  region.   The  flame  region 
was  located  within  the  developing  turbulent  boundary  layer  downstream  of 
the  reattachment  zone.   These  observations  are  in  qualitative  agreement  with 
experiment. 

A  comparison  of  predicted  plexiglas  fuel  regression  rates  with  experi- 
ment is  presented  in  Figure  9.   It  should  be  noted  that  the  location  of  the 
maximum  regression  rate  obtained  experimentally  was  affected  by  the  inlet 
air  distortion/turbulence  intensity.   The  model  predicts  a  maximum  regression 
rate  upstream  of  the  experimental  values  and  slightly  upstream  of  the  Krall- 
Sparrow"^  location  for  maximum  heat  transfer  in  non-reacting  flow.  The  aver- 
age fuel  regression  rate  was  in  good  agreement  with  experiment  and  the  regres- 
sion pattern  was  qualitatively  correct  (although  shifted  upstream) .   In  the 
present  model  only  convective  transfer  to  the  fuel  surface  is  considered.   The 
predicted  regression  profile  may  not  be  correct  for  fuel  systems  which  pro- 
duce large  amounts  of  radiative  heat  transfer  to  the  surface. 
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Figure  10  compares  the  predicted  centerline  turbulence  intensity 
(assuming  isotropic  turbulence)  with  experimental  data  for  non-reacting  flow 
The  behavior  with  axial  distance  was  generally  in  good  agreement  but  the 
inlet  region  exhibited  behavior  not  observed  experimentally.   The  model 
overestimates  the  velocity  increase  of  the  air  as  it  enters  the  combustor, 
which  in  turn  reduces  the  predicted  turbulence  intensity. 

Figure  11  compares  experimental  and  predicted  axial  pressure  profiles. 
The  experiment  and  theory  were  not  for  identical  conditions,  but  in  general 
it  is  seen  that  the  model  predicts  the  pressure  to  level-off  upstream  from 
experimental  measurements.   This  was  expected  since  the  model  underestimated 
the  distance  to  the  reattachment  point.   Pressure  calculations  were  very 
sensitive  to  solution  accuracy  since  they  depend  upon  second  derivatives  of 
stream-function. 

Figure  12  presents  predicted  radial  temperature  profiles  as  a  function  of 
inlet  step  size  and  air  flow  rate  at  a  location  downstream  of  flow  reattach- 
ment.  As  the  air  flow  rate  was  reduced  for  a  fixed  step  size,  the  peak 
temperature  moved  further  from  the  wall,  the  core  flow  became  hotter,  and 
the  overall  fuel-air  ratio  increased.   This  results  from  the  fact  that  regres- 
sion rate  varies  as  the  0.3  to  0.5  power  of  G  .   .   Actually  the  model  predicted 
that  f  ~  G  '      ,  whereas  experimentally  it  was  found  for  plexiglas  that  f  ~  G  ! 

The  recirculation  zone  in  the  solid  fuel  ramjet  acts  as  the  flame  sta- 
bilization device.   It  operates  fuel  rich  but,  when  operating  within  certain 
limits,  provides  enough  energy  to  heat  the  air  in  the  shear  layer  to  the 
point  where  the  downstream  boundary  layer  flame  can  be  stabilized  in 
the  high  velocity  flow.   As  the  inlet  flow  rate  was  increased  the  tempera- 
ture in  the  upstream  portion  of  the  recirculation  zone  was  found  to  increase. 
This  resulted  from  the  increased  air  flux  into  the  fuel  rich  portion  of  the 
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zone.  However,  this  same  increased  air  flux  into  the  recirculation  zone 
decreased  the  temperature  in  the  downstream  portion  of  the  zone  and  at  the 
reattachment  position.   Thus,  although  the  peak  temperatures  in  the  recircula- 
tion zone  increased  with  increased  inlet  air  velocity,  the  oxygen  reaching  the 
downstream  flame  zone  decreased  in  temperature.   In  addition,  the  boundary 
layer  flame  narrows  and  moves  closer  to  the  wall.   Eventually  a  condition  would 
be  reached  where  the  oxygen  reaching  the  boundary  layer  flame  region  would  be 
too  low  in  temperature  to  sustain  the  narrowed  flame  which  exists  in  the  rela- 
tively high  velocity  region  close  to  the  fuel  surface.   Earlier  work  with  the 

1  4 
model  '  ,  in  which  one-step  finite  rate  kinetics  where  employed,  also  demon- 
strated that  increased  inlet  velocity  can  cause  flame  blow-off  in  either  the 
recirculation  region  or  the  boundary  layer  region. 

Decreasing  the  inlet  step  size  decreased  the  length  of  the  shear  layer 
(reattachment  occurred  further  upstream)  and  decreased  the  air  entrainment 
into  the  recirculation  region.  Although  the  shear  layer  was  shorter  in 
length,  contact  time  with  the  recirculation  zone  did  not  necessarily  decrease 
because  the  inlet  velocity  was  lower.   However,  the  model  predictions  showed 
a  decreased  temperature  in  the  recirculation  zone,  resulting  from  a  still  more 
fuel  rich  mixture,  and  lower  temperatures  near  reattachment.  Again  it  would 
be  expected  that  a  flammability  limit  would  be  reached  with  increasing  inlet 
diameter. 

Figure  13  presents  the  calculated  percentage  of  unburned  fuel  as  a  function 

of  air  flow  rate.  As  discussed  above,  the  fuel  mass  fraction  was  the  least 

accurate  calculation.   Thus,  the  results  can  only  be  interpreted 

qualitative.   As   G  .    is  decreased,  an  increasing  amount  of  unreacted  and 
*  air 

hotter  fuel  is  passed  into  the  aft  mixing  region  downstream  from  the  end  of 
the  fuel  grain. 
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Schadow  has  made  measurements  of  temperature  and  species  concentrations 
in  a  reacting  solid  fuel  ramjet.  His  experiments  utilized  an  all-hydrocarbon 
fuel  which  results  in  significant  radiative  heat  transfer  to  the  fuel  surface 
and  higher  fuel  regression  rates „  However,  the  data  are  all  that  is  currently 
available  for  SFRJ  model  validation  in  the  reacting  flow  environment.  Because 
the  regression  rate  predicted  using  only  convective  heat  transfer  was  too  low 
for  the  radiative  system,  the  predicted  regression  rates  along  the  length  of 
the  fuel  grain  were  arbitrarily  increased  by  30%.  This  cannot  be  expected  to 
yield  the  correct  axial  variations  in  regression  rate  but  was  done  to  obtain 
a  more  realistic  value  of  the  blowing  parameter  for  the  higher  regression 
rate  system. 

Figs.  14,  15,  and  16  present  the  theoretical  temperature  profiles  compared 
to  the  data  of  Schadow  for  different  air  mass  fluxes  and  axial  locations 
within  the  fuel  grain.   The  theoretical  temperature  profiles  are  not  for 
exactly  the  same  Lp/Dp  locations  as  the  experimental  data.   The  profiles  at  the 
nearest  corresponding  grid  line  have  been  presented. 

Fig.  14  presents  data  for  the  high  mass  flux  condition.   The  model  pre- 
dicts more  complete  mixing  and  therefore  more  uniform  temperatures  in  the 
recirculation  zone  (Lp/D  =  0.52)  than  the  data.   In  the  region  near  reattach- 
ment (Lp/D  =  2.3)  the  experimental  data  indicates  an  apparent  "two-flame"  zone. 
It  is  possible  that  just  upstream  of  reattachment  one  flame  may  exist  very 
close  to  the  wall  and  another  could  result,  farther  from  the  wall,  from  the 
recirculating,  reacting  shear  layer.   However,  the  data  of  Schadow  has  not 
been  corrected  for  radiation  effects  and  an  L  /D  of  approximately  2.0  should 
be  the  most  downstream  location  of  reattachment  for  the  inlet  step  height 
employed.   If  this  "two-flame"  region  is  characteristic  of  the  flow  near 
reattachment,  then  the  model  fails  to  predict  this  detail.   There  could  be 
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several  reasons  for  this  possible  weakness  in  the  model.   Radiation  is 
neglected  in  a  system  which  produces  large  amounts  of  radiation.   In  addition 
the  K-£  model  and/or  the  boundary  conditions  specified  on  the  fuel  surface 
for  shear  stress  could  result  in  a  "smearing-out"  of  the  species  distribu- 
tions in  the  radial  direction  near  the  wall  where  large  property  gradients 
exist.  Downstream  of  reattachment  where  the  boundary  layer  develops 
(Lp/D  =  5.2,  6.1)  the  profiles  are  in  better  agreement  with  experiment.   The 
centerline  temperature  increase  with  axial  distance  is  reasonably  predicted. 

Similar  results  for  lower  air  fluxes  can  be  observed  in  Figs.  15  and  16. 
The  model,  without  radiation-convection  coupling  to  determine  the  fuel 
regression  rate,  predicted  a  stronger  dependence  of  regression  rate  on  air 
flux  than  obtained  experimentally.   This  will  also  affect  the  predicted 
temperature  profiles.   The  model  does  correctly  predict  the  increasing 
boundary  layer  thickness  with  resulting  peak  temperatures  further  from  the 
wall  and  higher  centerline  temperatures  as  the  air  flux  is  decreased.   The 
regression  rate  variation  with  axial  distance  is  different  for  the  all- 
hydrocarbon  system  than  for  the  convection  dominated  PMM  fuel  system.   The 
model  does  a  reasonable  job  of  predicting  this  profile  for  PMM.   For  both 
fuels  the  regression  rate  increases  from  the  head-end  to  near  the  reattachment 
zone  and  then  decreases.   Further  downstream  the  PMM  regression  rates  (and 
the  model  predictions)  continue  to  decrease  and/or  level  off  whereas  the 
all-hydrocarbon  regression  rates  begin  to  increase  again. 

The  model  assumed  mixing  limited  and  complete  combustion  so  that 
"products"  consisted  only  of  H_0  and  CO..   This  is  obviously  an  oversimplifi- 
cation but  perhaps  warranted  when  compared  to  the  other  simplifications 
employed  in  the  model. t     Schadow  has  measured  with  a  probe  the  radial 
variation  in  moles  of  free  CL  per  mole  of  N?  at  the  exit  plane  of  the  fuel 
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grain.   A  comparison  of  his  data  with  the  model  is  presented  in  Fig.  17. 
Schadow  found  that  in  order  to  predict  the  correct  variation  in  combustion 
temperature  with  mixture  ratio  in  the  fuel  port,  40%  of  the  carbon  present 
had  to  be  assumed  as  unreactive  soot  particles.   As  a  qualitative  attempt  to 
compare  the  model  (with  infinite  rate  kinetics  and  no  free  carbon)  with  the 
data  of  Schadow,  40%  of  the  carbon  in  "CO  products"  was  assumed  unreactive. 
The  oxygen  from  the  CO-  was  combined  with  the  unreacted  0-  to  produce  a  new 
profile  for  the  moles  of  "free"  0    .      The  results  are  presented  in  Fig.  17. 
and  show  reasonable  agreement  with  experiment. 

D.   Conclusions  and  Current  Work 

Temperature,  pressure,  turbulence  intensity,  and  species  profiles  as  well 
as  the  flow  reattachment  position  were  in  qualitative  agreement  with  experi- 
ment.  The  largest  discrepancies  existed  in  the  recirculation  region  and 
near  the  fuel  grain  at  reattachment.  Regression  rate  profiles  were  in  good 
agreement  except  for  the  all-hydrocarbon  fuel  far  downstream  from  reattachment. 
The  model  appears  to  be  adequate  for  qualitative  evaluation  of  the  effects  of 
geometry  and  test  conditions  on  the  flame  pattern  and  the  amount  of  unburned 
fuel  entering  the  aft  mixing  chamber.   In  order  to  obtain  better  quantitative 
results,  radiation  and  finite  rate,  multiple-step  kinetics  should  be  incor- 
porated into  the  model.   Further  work  is  required  to  validate  the  k-e  model, 
especially  in  the  near-wall  regions. 

The  Mf— oo  variables  make  it  difficult  to  obtain  accurate  predictions  of 
the  pressure  field  and  to  specify  boundary  conditions  on  the  walls.   The 
point  iterative  method  employed  also  requires  a  large  storage  space  and 
long  computation  times. 
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Current  work  with  the  model  includes  (a)  extending  its  application  to 
include  the  aft-mixing  geometry  as  shown  in  Fig.  1  and  (b)  using  the 
primary  variables  of  pressure  and  velocity  in  place  of  the  H'-co  variables 
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Fig.    1.      Schematic  of   Solid  Fuel  Ramjet 
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Fig.  2.   Schematic  of  Solid  Fuel  Ramjet  Apparatus 
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Fig.  3.   PMM  Regression  Rate-  Non- Bypass 
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Fig.    4.      PMM  Regression  Rate,    Bypass 
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Fig.  5.   PMM  Combustion  Efficiency  vs.  Air  Flux 
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Fig.  6.   Near-Wall  Control  Volumes 
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Fig.  8.   Predicted  Streamline  and  Temperature  Distributions 
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Fig.  11.   Axial  Pressure  Distributions 
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Fig.  13.   Predicted  Fraction  of  Unburned  Fuel  at  Grain  Exit 
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